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Recent results for the asymptotic behavior of fermion scattering amplitudes in the Sudakov limit are presented 
including next-to-leading logarithmic corrections. These are used for the analysis of the dominant electroweak 
corrections to the fermion-antifermion pair production in e + e~ annihilation at high energy. 



1. Introduction 

Four fermion processes are generally consid- 
ered as benchmark processes at high energy col- 
liders, with electron positron annihilation into 
muon or quark pairs at LEP and the Drell Yan 
process at hadron colliders as characteristic ex- 
amples. Within the presently accessible energy 
region, typically up to 200 GeV, radiative cor- 
rections are dominated by the shift in the W 
and Z masses as parametrized by the p parame- 
ter and by the running of the coupling constant. 
Vertex corrections and box diagrams involving 
gauge bosons are generally of minor importance. 
In the TeV region, accessible at future colliders 
like the LHC or TESLA, this picture changes 
drastically. A new class of effects starts to be- 
come relevant and rapidly dominant which are 
generally denoted as double logarithmic correc- 
tions and which were first observed by Sudakov 
[ |l| in the context of quantum electrodynamics 
for reactions with a tight cut on the radiated en- 
ergy of the photons. For electroweak interactions 
large negative corrections arise from the exchange 
of gauge bosons which remain uncompensated if 
one restricts the analysis to exclusive final states, 
consisting e. g. of a fermion antifermion pair only. 
The discussion of double logarithmic corrections 
is fairly straightforward for a theory with massive 
gauge bosons only. An important complication 
arises from the presence of massless photons in 
the final state. Events with soft and hard photon 
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radiation are normally included in the sample - 
whence a "semiinclusive" definition of the cross 
section is closest to the the actual experimental 
analysis. 

One loop corrections to the four fermion pro- 
cess are available since long (see e.g. [ || ) . In this 
case the separation of photonic and weak correc- 
tions is still possible and, employing axial gauge, 
the leading logarithms can be trivially attributed 
to fermion self energy diagrams with virtual W 
or Z boson emission [ || . 

Results for the leading logarithmic corrections 
in higher orders have been presented in [||, EL pU|] 
(for a more recent discussion see also [ M ISf). 
The first three results disagree, a consequence 
of different requirements on the inclusion or ex- 
clusion of hard photon emmission and a differ- 
ent treatment of the virtual corrections. In [ 
[| it was demonstrated that the approach de- 
veloped in [ |3j, where axial gauge was adopted 
and both real and virtual photon emmission ne- 
glected, is numerically nearly identical to the re- 
sults of [ |5| where virtual photon corrections and 
real photon emmission are included — a conse- 
quence of the smallness of the weak mixing an- 
gle, sin 2 6w — 0.23. Furthermore, an analysis of 
subleading logarithms for the form factor, for the 
four fermion scattering in a spontaneously bro- 
ken SU(2) gauge theory and, last not least, for 
the Standard Model with W, Z and the massless 
photon was performed in [ || . 

Before presenting a brief review of the tech- 
niques and some results from [ ^| let us empha- 
size the importance of subleading corrections. In 
the energy region relevant at LHC or TESLA, i.e. 
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around 1 to 2 TEV, this is already evident from 
the one loop vertex correction. In the timelike 
region 



\F B +6F\ 2 wF b (1 
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For characteristic values of the coupling, energy 
and mass of aw/v = 10~ 2 , Cp = 3/4 and 
s/M 2 = 10 2 one finds 



aw 

7T 



= 1.5 x 10~ 2 (-5.30 + 3.45-1.70) = 0.053. 



Large compensations between leading and sub- 
leading terms are observed and, in fact, this pat- 
tern will reappear for the Standard Model as dis- 
cussed below. 

The treatment of [ |(| is based on evolution 
equations that govern the dynamics of the am- 
plitudes in the Sudakov limit as obtained in [ 
H| [nj |ll|, |l2, 13 . In [ ^| this approach was applied 
to the next-to-leading analysis of the Abelian 
form factor and the four fermion amplitude in the 
SU(N) gauge theory. Functions that enter the 
evolution equations in the next-to-leading loga- 
rithmic approximation were evaluated by using, 
as an input, asymptotic expansions of one-loop 
diagrams. The solution of these equations lead 
to a summation of the leading and subleading 
Sudakov logarithms. The expansion of one-loop 
diagrams through the so-called generalized strat- 
egy of regions [ [l4| (see also [ |l5)) identifies in 
a systematic way the nature of various contribu- 
tions and the origin of logarithms. This strategy 
is based on expanding integrands of Feynman in- 
tegrals in typical regions and extending the inte- 
gration domains to the whole space of the loop 
momenta so that a crucial difference with respect 
to the standard approach [ ||, 10, [0], |l2|, [l3j is 
the absence of cut-offs that specify the regions in 
individual terms of the expansions. 

2. The Abelian form factor in the Sudakov 
limit 

Let us first analyse the (vector) form factor 
which determines the amplitude of the fermion 



scattering in the external Abelian field. In Born 
approximation 



We consider the limit s = (pi — P2) 2 



(2) 



with on-shcll massless fermions, p\ = p\ — 0, and 
gauge bosons with a small non-zero mass M 2 <C 
— s. For convenience pi ; 2 = (Q/2, 0, 0, =fQ/2) so 
that 2pip2 = Q 2 = — s. 

The asymptotic behaviour can be found by 
solving the corresponding evolution equation [ 0] 
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For the non-Abelian gauge theory, this equation 
was first derived in [ [l0| by factorizing collinear 
logarithms in the axial gauge. Its solution is 



^ = F (a(M 2 ))exp^ / — 
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For a proper treatment of the next-to-leading log- 
arithms one must keep renormalization group cor- 
rections to the leading logarithmic approxima- 
tion as well as single infrared and renormalization 
group logarithms. In this approximation 



F = F (a) exp 



^ d i f ^r^ xl)) (5) 

M 2 x JM 2 x 
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The leading terms of the functions 7, £ and £ are 
obtained from the one loop analysis and the one 
loop running of a in the argument of the function 
7 should be taken into account. 

In the covariant gauge, the self energy inser- 
tions to the external fermion lines do not give 
Q-dependent contributions. The one loop calcu- 
lation of the vertex correction gives 

F = ^-C F (-Vo + 21/i + 2(1 - 2e)V 2 ~ V±) F B , 

(6) 
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where Cp is the quadratic Casimir operator of the 
fundamental representation and the functions Vi 
are obtained from 



d d k 
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To expand these integrals in the limit Q 2 ^> 
M 2 we apply a generalized strategy of regions for- 
mulated in [ and discussed using characteristic 
two- loop examples in [ [l5| : 

• Consider various regions of the loop mo- 
menta and expand, in every region, the in- 
tegrand in Taylor series with respect to the 
parameters that are there considered small; 

• Integrate the expanded integrand over the 
whole integration domain of the loop mo- 
menta; 

• Put to zero any scaleless integral. 

The following "typical" regions arise in the Su- 
dakov limit [ ll(| : 



hard (h): 

1- collinear (lc): 

2- collinear (2c): 

soft (s): 
ultrasoft (us): 



Q 

M 2 /Q, k ll2 



M 



k~Q 
fc + k 3 
k - k 3 
k - k 3 ~ Q 

fc + fc 3 -M 2 /Q, fci, 2 ~M 

fe~ M 

k ~ M 2 /Q . 



Keeping the leading power in the expansion in 
the limit Q 2 /M 2 — > oo one observes [ |) that the 
leading double logarithm results from the hard 



region, whereas the single logarithm receives con- 
tributions from the hard and collinear regions as 
well. Soft regions do not contribute, at least in 
the leading power. Combining the vertex correc- 
tion with the Q 2 independent fermion self ener- 
gies one arrives at a finite result. 

From the one-loop result one derives 



7(a) = -Ci 



a 
2^ 
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The total double logarithms originate from the 
hard region. This determines the scale of the 
coupling constant in the second order logarith- 
mic derivative of the form factor in Q. At the 
same time we cannot distinguish, in the one loop 
approximation, the contribution to the functions 
C and £ coming from the collinear region because 
this region includes both Q and M scales. For 
the sum of these functions we find 



(9) 



Finally, in the NLO logarithmic approximation 
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in agreement with the result of [ [lj, Q . 



3. The four fermion amplitude 

We study the limit of the fixed-angle scatter- 
ing when all the invariant energy and momen- 
tum transfers of the process are much larger 
than the typical mass scale of internal particles 
|s| ~ \t\ ~ \u\ 3> M 2 . Besides the extra kine- 
matical variable the analysis of the four fermion 
amplitude is more complicated by the presence 
of different "color" and Lorentz structures. The 
Born amplitude, for example, can be expanded in 
the basis of color/chiral amplitudes 

A B = !l A ^ = ^T F ^±(At L + At R ) (11) 
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where 
A x 
A d LL 

A LR 



i>2 (P2 )t a l^l (pl ) ^4 (P4)*°7 M V>3 (?3 ) » 

^alT/^ii^iT/^i > ( 12 ) 



and so on. Here t a is the SU(N) generator, p\, 
P3 are incoming and P2, Pa outgoing momenta so 
that t — (pi ~Pi) 2 and u = (p% +P3) 2 = —(s + t). 
For the moment we consider a parity conserving 
theory. Hence only two chiral amplitudes are in- 
dependent, for example, LL and LR. Similarly 
only two color amplitudes are independent, for 
example, A x and A d . 

Let us first compute the one loop corrections, 
employing again the strategy of regions. The to- 
tal contribution from vertex and box diagrams 
in the logarithmic approximation is independent 
from chirality and the same both for the LL and 
LR amplitudes: 
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and the same both for the LL and LR amplitudes. 

Now the collinear logarithms can be separated 
from the total one-loop correction. For each 
fermion-antifermion pair, they form the exponen- 
tial factor found in the previous section (eq. (|Tc|) ) . 
This factor in addition incorporates the renormal- 
ization group logarithms which are not absorbed 
by changing the normalization scale of the gauge 
coupling. The rest of the single logarithms in 
eq. ( |l3| ) is of the soft nature. Let us denote by A 
the amplitude with the collinear logarithms fac- 
tored out. It can be represented as a vector in the 
basis A x , A d and satisfies the following evolution 
equation [ |ll|, [r| 

^i = xMQ 2 ))i, (14) 



where x 1S the matrix of the "soft" anomalous 
dimensions. From eq. ( |l4| ) we find the elements 
of this matrix to be, in units of a/(4ir), 
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The solution of eq. (fll) reads 



I = A°(a(M 2 ))exp 



J 4°(a(M 2 ))exp 



da; 



M 2 x 

M 2 X 



Xi(a(x)) 



(16) 



X2{a(x)) 



where Xi are eigenvalues of the x matrix and A® 
are Q-independent vectors. 

/From the asymptotic expansion of the box dia- 
grams one finds that only the hard parts con- 
tribute to eq. (|l4|). This fixes the scale of a in 
this equation to be Q. 

In the Abelian case, there are no different color 
amplitudes and there is only one anomalous di- 
mension 

X = 41n(|). (17) 

4. Sudakov logarithms in electroweak pro- 
cesses 

We are interested in the process /'/' — > //. In 
the Born approximation, its amplitude is of the 
following form [ pf 



A, 



where 



I, J=L, R V - 



A IJ f = fllvf'lfjlufj . 



(18) 
(19) 



To analyze the electroweak correction to the 
above process we use the approximation with the 
W and Z bosons of the same mass M and mass- 
less quarks and leptons. The photon is massless, 
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and the corresponding infrared divergent contri- 
butions should be accompanied by the real soft 
photon radiation integrated to some resolution 
energy u> res to get an infrared safe cross section 
independent on an auxiliary photon mass. At the 
same time the massive gauge bosons are supposed 
to be detected as separate particles. In practice, 
the resolution energy is much less than the W 
(Z) boson mass so the soft photon emission is 
of the QED nature. This cancels the infrared 
singularities of the QED virtual correction. We 
should therefore separate the QED virtual cor- 
rection from the complete result computed with 
the photon of some mass A and then evaluate the 
QED virtual corrections together with the real 
soft photon radiation effects with vanishing A. It 
is convenient to subtract the QED contribution 
computed with the photon of the mass M from 
the obtained result for the virtual corrections and 
then take the limit A — > for the sum of QED 
virtual and real photon contributions to the total 
amplitude. In the language of the approach of 
ref. [ H , this prescription means that we use the 
auxiliary photon mass A as a variable of the evo- 
lution equation below the scale M and the sub- 
traction fixes a relevant initial condition for this 
differential equation. This leads to a modification 
of the collinear factor and the soft anomalous di- 
mensions. 

Thus we keep always a cutoff ui res <C M and 
do not display the QED Sudakov factor arising 
from this suppression of photon radiation. The 
remaining "electroweak" universal collinear factor 
for each fermion-antifcrmion pair becomes 



exp 



T f (T f + l)+t 2 w -L- s 2 w Q) 



x (L(s) -31(a)) 



(20) 



The soft anomalous dimension for / and/or J = 
R is Abelian and, in units of <? 2 /(167r 2 ), reads 

X = {t 2 w Y f ,Y f - 4a 2 w Q f ,Q f ) In (|) , (21) 

and the matrix of the soft anomalous dimension 
for / = J = L is 



XAA = 



XdX 
Xdd 



+ {twYf'Yf - 4s 2 w Q r Q f + 2) In (|) 

!>»(!) 



= 4 In 
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(twYf'Y f -4s 2 w Q f ,Q f ) ln(- 



The one-loop leading and subleading logarithms 
can be directly obtained from eq. (|l3|). The two- 
loop leading (infrared) logarithms are determined 
by the second order term of the expansion of 
the double (soft x collinear) logarithmic part of 
the collinear factors (20). The two-loop next- 
to-leading logarithms are generated by the inter- 
ference between the first order terms of the ex- 
pansion of the double (soft x collinear) and sin- 
gle (soft+collinear+renormalization group) log- 
arithmic exponents and can also be found in [ 
§. With the expression for the chiral amplitudes 
at hand, we can compute the leading and sub- 
leading logarithmic corrections to the basic ob- 
servables for e + e~ — * //. Let us, for example, 
consider the total cross sections of the quark- 
antiquark//i + /i~ production in the e + e~ annihi- 
lation. In the two loop approximation, the lead- 
ing and next-to-leading Sudakov corrections to 
the cross sections read 

a/a B (e + e- -> QQ) = 1+ 5.301(a) - 1.66£(s) 

-12Ml(a)L(a) + 1.92L 2 (s) , 
<j/<j B (e + e- -*qq)=l + 20.541(a) - 2.Y7L(a) 

-53.72l(s)L(s) + 2.79L 2 (s) , 

cr/cr B {e + e- -> = l+10.09Z(s)-1.39L(s) 

-21.66Z(s)£(s) + 1.41£ 2 (s) , 

(23) 

where Q = u,c,t, q = d, s, b. Numerically, L(a) = 
0.07 (0.11) and 1(a) = 0.014 (0.017) respectively 
for y/a = 1 TeV and 2 TeV. 

Clearly, for energies at 1 and 2 TeV the two 
loop corrections are huge and amount up respec- 
tively to 5% and 7%. There is a cancellation 



6 



between the leading and subleading logarithms 
and for the above energy interval the sublead- 
ing contribution even exceeds the leading one. 
The higher order leading and next-to-leading cor- 
rections however do not exceed 1% level. They 
can be in principle resummed using the formu- 
lae given above. The leading and subleading cor- 
rections to the left right and forward backward 
asymmetries [ || are typically smaller. 

Our result for the one loop double logarithmic 
contribution is in agreement with [ . However 
the result for the one loop single infrared logarith- 
mic contribution differs from [ [l8| . The reason is 
that, in [ ^8), only the diagrams with heavy vir- 
tual bosons have been taken into account. There 
is an infrared safe contribution of the diagram 
with the virtual massless photon where the heavy 
boson mass serves as an infrared regulator that 
should be taken into account to get a complete 
(exponential) result. In one-loop approximation, 
this contribution comes from the box diagrams 
with the photon and Z boson running inside the 
loop[|l|. 

Our result for the two-loop double logarithmic 
contribution is in agreement with [ Js). On the 
other hand, the coefficients in front of the two- 
loop leading logarithms in eq. (|2^) with a few 
percent accuracy coincide with the result of [ |j| 
where the photon contributions were not consid- 
ered. This is related to the fact that the virtual 
photon contribution not included to the result of 
[ ||] is suppressed by a small factor s^- 
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